Astronomy & Astrophysics manuscript no. 13412 


©ESO 2010 


September 29, 2010 





The self-cohering tied-array 



P. Fridman 



ASTRON, Oudehoogeveensedijk 4, Dwingeloo, 7991PD, The Netherlands 
e-mail: fridman@astron.nl 



ABSTRACT 



o 

< 



Context. Large radio astronomy multi-element interferometers are frequently used as single dishes in a tied-array mode when signals 
from separate antennas are added. Phase shifts arising during wave propagation through a turbulent atmosphere can significantly 
reduce the effective area of an equivalent single dish. 

Aims. I aim to give estimates of the impact of the ionosphere and troposphere on the effectiveness of a radio interferometer working 
in tied-array mode. 

Methods. Statistical estimates of the effective area are calculated and the power-law of turbulent atmosphere irregularities has been 
used. A simple method of tied-array calibration using optimization techniques is proposed. 

Results. The impact of phase errors on the effectiveness of tied-arrays are given for low and high frequencies. Computer simulations 
demonstrate the efficacy of the proposed calibration algorithm. 
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1. Introduction 

Large radio astronomy multi-element interferometers (VLA, 
WSRT) are frequently used in the tied-array mode where sig- 
nals from separate antennas are added (Thompson et al. 2001, 
ch. 9.9). The output sum signal can be used in VLBI, pulsar and 
transients observations (Cordes 2004), SETI signals detection 
and the direct-to-Earth (DtE) reception of signals from cosmic 
apparatus (Jones 2004). In all these cases a radio interferometer 
works as a single-dish antenna with one output. Partial signals 
from antennas are properly phased to collect emission from a 
point-like radio source in the sky and track it during its siderial 
movement. Standard calibration procedure using a correlator is 
employed to provide the necessary phase corrections for each 
individual antenna. Random phase perturbations such as phase 
shifts arising during wave propagation through the turbulent at- 
mosphere can occur in the course of such observations. These 
phase errors reduce the total effective area of the tied-array and 
must be compensated for in real time. Although it is possible to 
store baseband data for processing after observations have taken 
place, the amount of data to be stored places a firm limit to the 
number of antennas that can be used in this manner. 
New large scale projects such as SKA and LOFAR will also be 
operating in tied-array mode. The impact of ionospheric and tro- 
pospheric phase errors on the tied-array is calculated in this pa- 
per. A simple method of correcting these errors using the output 
signal of the tied-array is also proposed here. 

2. Tied-array with random phase errors 

Voltage produced by the source at the output of the planar tied- 
array is 



r Sin -is the distance between the source and n-th array element, 
n a is the number of antennas in the array, 5 n is the instrumental 
phase shift. The distance r Stn can be represented as the module 
of the vector difference s - p w , where p„ is the position vector of 
n-th array element: 

r s ,n = ^\S~Pn\ 2 = ((S x - p n , x f + (s y - p n , y f + ^) 1/2 = 

(r 2 - 2p n , x s x - 2p n , y Sy + pl x + ply) 111 *r s - 



Pn,x^x Pn,ySy 



r s ~ Pn,x cos(or JfJC ) - p n , y cos(q^) = r s - p„u s {2) 

s X9 s y ,s z are the components of vector s and p n , x ,Pn, y are the 
components of vector p w , u s = [cos(o^, x ), cos(a Sr y)]. 
With these new notations (1) can be rewritten as: 

n a 

E(u s ) = exp(j2nft - jkr s ) ^ a^ n j\_kQ n u s + 6 n ]). (3) 



An array is directed at s when these conditions are satisfied: 

- kp n u s +6 n = 0,n= \...n a . (4) 

The manifold of signals received in other directions (u ^ u s ) 
forms the field pattern: 



£(u, u s ) = Qxp[j27rft - jk(r s - r u )] ^ a^ n exp[-j£p w (u, - u)]. 

n=\ 

(5) 

The power pattern is the expected value of the product 



P(u,u s ) =< E(u, u s )E(u, u s ) >= 



n a n a 



E(s) = expOWf) V a s , n exp[-j(kr s , n + 6 n )] (1) < h 2>^»- exp[-jfe(p„ - Pm )(« f - u)]} > 

' * 77=1 m— 1 



(6) 



where s is the source vector, k = 2n/A, A is the wavelength, / is 
the frequency, a s , n - is the signal amplitude at n-th array element, 



where () denotes a complex conjugate and <> denotes a time av- 
erage. 
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In the absence of phase errors and in the direction of the radio 
source, i.e., for u = u s ,a u , n = a,P(u s ,u s ) = n 2 a a 2 . The dc com- 
ponent proportional to the sum of the system noise power at each 
antenna is omitted here and considered to be a constant value and 
therefore not relevant. 

In the presence of phase errors produced by, for example, the at- 
mosphere, there are additional phase terms 6 n ,atm ^ in (6) and 
the power received in the direction is: 

Ps am (u s ,u s ) = 

n a n a 

^ X X < eX PW<W™ ~ $m,atm) > = 
n- 1 in- 1 

a ZjZj exp[ 2 — ] ( } 

n- 1 m- 1 

where D^ atm {b mn ) is the variance of the random phase difference 
for the baseline b mn = |p w - p m |: 

Ds,atm(bmn) =< (<>n,atm ~ <>m,atm) 2 > (8) 

It is assumed that the random phase difference A = 6 nAtm - S nuitm 
has normal probability distribution with zero mean and variance 
<x 2 . In this case (7) was obtained using this relation: 

1 r°° 

< exp(jA) >= exp(jA)exp(-A/2o- 2 )dA = 

exp(-<7i/2) (9) 
Therefore, the loss produced by the phase errors is equal to: 

The variance D^ atm (b mn ) is the structure phase function of the 
turbulent atmosphere. The power-law (Kolmogorov) model will 
be used in the following sections to describe D^ atm (b mn ) both 
for the ionosphere and troposphere phase fluctuations (Tatarskii 
1978). 

2.1. Ionosphere 

Phase fluctuations due to the irregular spatial distribution of 
the refraction index during wave propagation through the iono- 
sphere are described with the power-law model of the turbulence 
spectrum. The electron density N in the ionosphere, considered 
as a function of spatial coordinates, has variations which are 
characterized by a structure function of electron density D N (b) 
(Thompson et al. 2001, ch. 13): 

D N (b) = C 2 N b\ (11) 

where y = 2/3, C 2 N is the constant, b is the baseline. This formula 
can be rewritten for the structure function of the refraction index: 

D n (b) = C 2 n b\ (12) 

where C 2 = ^rCfj, r e = 2.82 • 10" 15 m (electron radius), A is the 
wavelength. Finally, the ionosphere phase structure function is: 

D ion (b) = 2.91/c 2 C 2 n hb 5/3 = 2.91r 2 A 2 C 2 N hb 5/ \ (13) 

where h is the total propagation length through the irregularities 
of the ionosphere, A, h and b must be substituted in meters. The 
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Fig. 1. a) upper panel: square root of the ionosphere structure 
function of electrical length, in cm; N = 10 12 m" 3 , A(N)/N = 
0.01 and Lq = 10km calculated for three frequencies: 50MHz, 
100MHz and 200MHz; b) lower panel: Fried length for different 
values ofN,Lo,h. 

value C 2 N can be estimated from (11) assuming that the iono- 
sphere irregularities of electron density A(AT) have a maximum 
dimension equal to Lq\ 

A(N) 2 = C 2 N Lf. (14) 

For example, 

for A(AT)/W = 0.01 and N = 10 12 m" 3 (day time) we have for 
L = 10£m, C 2 = 2.154 • 10 17 m- 20/3 , 

for A(N)/N = 0.03 and L = 30/cm,C 2 = 9.322 • 10 17 m- 20/3 . 
Figure la shows the square root of the ionosphere structure func- 
tion for N = 10 12 m" 3 , A(N)/N = 0.01 and L = 10km, h = 
300km calculated for three frequencies: 50MHz, 100MHz and 
200MHz. 

Phase fluctuations can be also characterized by the Fried length: 
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r Fried = 3. 18 Jo, 



(15) 



where do is the baseline at which ^D ion = Ira d. For the parame- 
ters used in Fig. la r Fried = 1.16km for 50MHz, r Fried = 2.66km 
for 100MHz and r Fried = 6.1 Urn for 200MHz. Fig. lb shows 
how Fried length depends on the frequency for the different val- 
ues ofN,Lo,h. 

The minimal time interval at which it is necessary to re- 
peat calibrations can be calculated from Fried length: t ca \ - 
* 'Fried /(3.18v wO where v w ^ is the wind velocity. Thus, for ex- 
ample, for rpried = 6.11km and v w i n d - 50m/ sec, t ca \ - 38sec. 

Now the loss produced by ionosphere random phase errors 
can be calculated in the example of the array whose configura- 
tion is shown in Fig. 2a. It is the random planar 100-element 
array with the coordinates %i and yi represented by random nor- 
mal values with zero mean and standard deviation S C. For the 
array shown in Fig. 2a, SC = 1000m, therefore the maximum 
baseline is « 5000m. The distribution of baselines (histogram) 
is shown in Fig. 2b . Phase errors are maximal for the largest 
baselines but their relative number is less than the medium size 
baselines, therefore the signal loss for the array should take this 
particular distribution of baselines into account. Fig. 3 demon- 
strates the dependance of loss versus array size 5 SC. The curves 
are calculated for three frequencies: 50, 100 and 200MHz. 

2.2. Troposphere 

Phase fluctuations due to the irregular spatial distribution of 
the refraction index during wave propagation through the tropo- 
sphere are also described in the frame of the power-law spectrum 
turbulence model. The troposphere phase structure function is 
(Stotskii 1973, Carilli et al. 1999): 



D trop (b) = 2.91k 2 C 2 b 5/3 ,L <b<L l 
= 2.91& 2 C?Z? 2/3 ,Li <b<L 2 



2.91k 2 C L ,L 2 <b, 



(16) 



where Lo and L\ are the internal and external scales, respec- 
tively, of the isotropic three-dimensional turbulence model, Lo = 
0.1 - 1cm, U = 5.6Km and L 2 = 2000 - 3000/^m, the latter is 
determined by global meteorological variations. 
Factors C 2 and C 2 depend on the local content of water va- 
por and oxygen in the troposphere (weather conditions) and the 
values chosen for the purpose of calculation are C 2 = 6.23 • 
10- n m 1/3 and C 2 = 3.64 • 10" 7 m 4/3 . 

Fig. 4a represents the structure function (16) and Fig. 4b shows 
the Fried length as a function of the frequency. Fig. 5 demon- 
strates the dependance of loss versus array size 5 -S C. The curves 
are calculated for three frequencies: 1400, 5000 and 8400MHz. 

3. Self-cohering 

Observations in the tied-array mode (VLBI, transients, DtE) 
are pre-planned at any time and it is impossible to postpone 
them in order to choose better atmospheric conditions (for ex- 
ample, night time in the case of the ionosphere). The effective- 
ness of the synthesized aperture must be maximal during obser- 
vations which means that periodical calibrations are necessary. 
Traditional methods of radio interferometer calibration can be 
applied using the grid of calibration point sources. This calibra- 
tion must be made in real time with the help of available correla- 
tors which must work in parallel with the tied-array adder. Here, 
another method is proposed which uses the points in the direct 
images of the field-of-view with calibration sources. 
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Fig. 2. a) upper panel: random array configuration; b) lower 
panel: histogram of the baselines. 

It is presumed that a full calibration with the correlator 
has already been performed before each tied-array observation. 
During subsequent observations the total power output of the 
tied-array is used as a tracking tool and the proposed algo- 
rithm will introduce small phase corrections at short time inter- 
vals therefore keeping the amplitude of the calibration source 
at a prescribed level. This has similarities to the approach of 
(Muller&Buffington 1974) which is also described in (Tyson 
1991). 

The choice of calibration sources is the same as in traditional 
methods. 

Equation (7) corresponds to the synthesized beam when 
there are phase errors 5 n ,atm produced during propagation 
through the turbulent atmosphere. To eliminate 6 n , atm , compen- 
sation phase shifts d n ,com P are introduced at each ft-th array ele- 
ment. The phase of the signal corresponding to the direction 
at the n-th array element is 



6 n (u s ) = (2n/A)(-p n u s ) + 6 n + 6 nAtm - d n + 



(17) 
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Fig. 3. Loss produced by phase fluctuations in the ionosphere 
calculated for three frequencies: 50, 100 and 200MHz. The 
structure function from Fig. la is used. 

where d n ,com P is the compensation phase shift introduced for the 
correction of S n , atm - atmospheric phase error. The value of the 
signal power in the prescribed direction u s in the image with 
one or several calibration sources available in the field-of-view 
(FoV) can be obtained by convolution of the sky intensity Bq(u) 
with the synthesized beam P , com ^(u J ) 

BcompiUs) = B (U) Pcompfrs), (18) 

where © denotes convolution. The similar approach to imag- 
ing procedure with an array is considered in the direct imaging 
scheme, (Wright 2004). Having the output of the total power de- 
tector B comp (u s ) for the direction we can maximize this value 
by varying S n , comp , i.e., we have an optimization problem for n a 
values of compensation phases. In general, for each direction u s 
this problem can be written as 

8njomp(u s ) = arg max[B comp (u s )] 9 n = \...n a (19) 

3n,comp 

This scheme is illustrated in the following example of com- 
puter simulation. The 60-element spatially random planar array 
is represented in Fig. 6. The phase errors introduced in each of 
the array element signals are modeled by the two-dimensional 
random value, Fig. 7, with circular symmetrical spatial spec- 
tral density which decreases radially according to the (-11/3) 
power law. The instantaneous sample of phase errors at fre- 
quency 100MHz as a function of baseline length is given in Fig. 
8. This high level of phase errors has been chosen to demonstrate 
the effectiveness of the proposed method. 
The image containing three point sources is represented in Fig. 
9 (left panel) and the synthesized image in the presence of the 
phase errors (Fig. 8) is shown in Fig. 9 (middle panel), (iso- 
planicity being presumed). 

The value of the synthesized image in the direction of the 
largest source (lower left in the image) was used as the cost 
function. The genetic algorithm(GA) was applied because of 
the strong multi-modality of the cost function (19) and this al- 
gorithm finds the global maximum successfully. Genetic algo- 
rithms search for a solution to a set of variables by the use 
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Fig. 4. a) upper panel: Square root of the troposphere structure 
function of electrical length, in cm; b) lower panel: Fried length 
as a function of frequency. 



of simulated evolution, i.e., the survival of the fittest strategy. 
In contrast to calculus-based algorithms (conjugate gradients 
and quasi-Newtonian methods), GA, first introduced in (Holland 
1975), exploit a guided random technique during optimization 
procedure (Goldberg 1989, Michalewicz 1994, Charbonneau 
1995). 

GA optimizers are particularly effective when the goal is to find 
an approximate global maximum in a high-dimension, multi- 
modal function domain in a near-optimal manner. They are also 
largely independent of the starting point or initial guess. There is 
parallelism which allows for the exploitation of several areas of 
the solution space at the same time. This parallelism can be very 
useful in the implementation of GA on the multi-core platform 
and FPGA. In this article, computer simulation has been done on 
a PC (Intel Pentium, 2.5 Ghz, 1GB RAM) using Matlab 7.6.0. 
Specific GA operations (selection, crossover and mutation) have 
taken approximately 3% of the total computing time: 150 sec 
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Fig. 5. Loss produced by phase fluctuations in the troposphere 
calculated for three frequencies: 1400, 5000 and 8400MHz. The 
structure function from Fig. 4a is used. 

for 100 iterations (each iteration is the full cycle of these opera- 
tions). The rest of the computing time was spent on the calcula- 
tion of the cost function (formation of the beam with corrected 
phases - > convolution with the image - > total power output). 
But these calculations are necessary only in computer simula- 
tions: in reality the values of the cost function are supplied by 
the tied-array itself ("Nature" does the job). 

After applying the optimization procedure and introducing 
the resulting compensation phases the corrected image is shown 
in Fig. 9 (right panel). 

The contour presentations in Fig. 10 correspond to the undis- 
torted image (left panel), the image with phase errors (middle 
panel) and the image after correction (right panel), respectively. 
The corresponding synthesized beam is restored up to 0.94 of its 
undistorted value. 

There are some peculiarities in image processing with non- 
planar arrays (Perley 1999) but the tied-array mode concerns 
only point-like sources. Therefore, there is no difference in pla- 
nar and non-planar arrays in the context of this article (phase 
irregularities due to atmospheric turbulence), especially for the 
adaptive calibration procedure described here. 

4. Conclusions 

1 . The effective area of tied arrays may be significantly reduced 
by ionospheric and tropospheric phase irregularities at low 
and high frequencies, respectively. 

2. Observations are made at times (VLBI, transients monitor- 
ing, DtE) when it is impossible to choose quiet atmospheric 
conditions and real-time calibration is necessary and has to 
be fulfilled in parallel with observations. 

3. The total power at the auxiliary outputs of the tied-array, 
phased in the direction of calibration sources, can be used on 
a level with traditional calibration methods using correlators. 
Multi-beam facilities are necessary for creating these aux- 
iliary outputs. Optimization algorithms (genetic algorithms, 
simulated annealing) can be used to compensate for propa- 
gation phase errors by maximizing the amplitude of a chosen 
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Fig. 7. Spatial phase error distribution, projected on the array 
plane. 
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Fig. 8. Phase errors as function of baseline length. 

calibration source. The tied array can preserve its correctly 
phased state during lengthy observations using one or sev- 
eral auxiliary outputs, therefore working in the self-cohering 
regime. The proposed scheme does not exclude traditional 
methods of calibration - it is complementary to them. 

Acknowledgements. I am grateful to Roy Smits whose comments were very 
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image, no phase errors 



image, with phase errors 



image after correction 






Fig. 9. Left panel: synthesized image without phase errors; middle panel: synthesized image with phase errors; right panel: synthe- 
sized image after correction. 
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Fig. 10. Contour presentations, left panel: synthesized image without phase errors; middle panel: synthesized image with phase 
errors; right panel: synthesized image after correction. 
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